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ABSTRACT 

We have improved upon the method of smoothing supernovae data to reconstruct the 
expansion history of the universe, h(z), using two latest datasets, Gold and SNLS. 
The reconstruction process does not employ any parameterisation and is independent 
of any dark energy model. The reconstructed h(z) is used to derive the distance factor 
A up to redshift 0.35 and the results are compared with the given value of A from 
detection of baryon acoustic oscillation peak (BAO). We find very good agreement be- 
tween supernovae observations and the results from BAO for f2o m w 0.276±0.023. The 
estimated values of fiom are completely model-independent and are only based on ob- 
servational data. The derived values of f^om are then used to reconstruct the equation 
of state of dark energy, w(z). Using our smoothing method we can demonstrate that 
while SNLS data are in very good agreement with ACDM, the Gold sample slightly 
prefers evolving dark energy. We also show that proper estimation of the equation of 
state of dark energy at the high redshifts would be impossible at the current status of 
observations. 
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1 INTRODUCTION 

Over the last two decades, cosmology has entered the stage 
of 'precision' science, which involve significant improve- 
ments in observational techniques, implementation of more 
powerful statistical and mathematical tools, and, of course, 
greatly advanced computational facilities. While these ad- 
vancements have yielded much new insight into the subject, 
many important questions still remain unanswered. The na- 
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Copeland fc Sami fc Tsuiikawal 2004 ). Supernovae data, 
which gave the first indication of the accelerated expansion 
of the universe, are expected to elucidate this interesting 
question further, as the quality of the data steadily improves 
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Fay fc Tavakoll 12006 ). A comprehensiv e rece nt review has 
been gi ven by ISahni fc Starobinskvl |2006). An earlier 
work by Shafie loo et all (2006) suggested a non-parametric 
method based on smoothing the supernova data over 
redshift in order to reconstruct cosmological quantities, 
including the expansion rate, h(z), and the equation of 
state of dark energy, w(z), in a model-independent manner. 
In this approach, the data are dealt with directly, and one 
does not rely on a parametric functional form for fitting any 
of the quantities dL(z), h(z) or w(z). The result obtained 
by using this approach is, therefore, expected to be model- 
independent. This method was shown to be successful in 
discriminating between different models of dark energy if 
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the quality of data is commensurate with that expected 
from the future SuperNova Acceleration Probe (SNAP). In 
this paper we improve the smoothing method and apply it 
to tw o recent set s of supernovae d ata: Gold IjRiess et al.l 
l2006h and SNLS (|Astier et al.ll2005l l. We then compare the 
derived expansion history of unive rse with the results o f 
baryon acoustic peak observations l|Eisenstein et al]|2005h . 
Specifically, we use the improved smoothing method to 
reconstruct the Hubble parameter, h(z), and then derive 
the distance factor, A, up to a redshift of 0.35 independently 
of the assumption of any cosmological model. This derived 
value, based on supernovae data, is then compared with 
the distance factor A (which is also claimed to be relatively 
independent of dark energy model) being determined by 
the detection of the baryon acoustic oscillation peak. One 
of the main results of this paper is that there is a good 
agreement between supernovae data (both Gold and SNLS) 
and baryon acoustic peak observations for the values 
fiom ~ 0.276 ± 0.023. The derived value of fiom is then 
used to reconstruct the equation of state of dark energy for 
both supernovae datasets. We should emphasize here that 
all the results in this paper are only based on observational 
data and no theoretical model has been assumed. This is 
an advantage of this method over the functional fitting 
methods in which the reconstructed results are biased by an 
assumed functional form or a theoretical model. The paper 
is organised as follows. In Section II we briefly explain the 
smoothing method and we estimate the accuracy of the 
method based on the quality and the quantity of current 
datasets. In Section III we apply the smoothing method on 
the Gold supernovae dataset and by using the results of 
detection of baryon acoustic oscillation peak, we reconstruct 
w[z). In Section IV, we carry out a similar treatment on the 
SNLS dataset. Finally, in Section V we discuss our results 
along with some concluding remarks. 



2 METHOD OF SMOOTHING 

The method of smoothing belongs to the category of non- 
parametric methods of reconstruction which is complemen- 
tary to the approach of fitting a parametric ansatz to the 
dark energy density or the equation of state. Most pa- 
pers using the non-parametric approach exploit a kind of 
top-hat smoothing in redshift space. Instead, we follow 
a procedure which is well known and frequently used in 
the analysis of large-sc ale structure |Coles fc Lucchinll 19951 ; 
iMartinez k, Saarll2002h ; namely, we attempt to smooth noisy 
data directly using a Gaussian smoothing function. In this 
method we apply Gaussian smoothing to supernova data 
(which are of the form {\nd,L(zi), z,}) in order to extract 
information about important cosmological parameters such 
as H(z) and w(z). The smoothing algorithm calculates the 
luminosity distance at any arbitrary redshift z to be 



lnd L (z,A) s = ln d L (z) 9 +iV(2)^[mdx,(zi) -In d L ( Zi 
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Here lndL(z,A) s is the smoothed luminosity distance 
at any redshift z which depends on luminosity distances of 
each SNe event with the redshift Zi, and N(z) is a nor- 
malisation parameter. The quantity \ndL(z) 9 represents a 
guess background model which we subtract from the data 
before smoothing. This approach allows us to smooth noise 
only, and not the luminosity distance. After noise smooth- 
ing, we add back the guess model to recover the luminosity 
distance. This procedure is helpful in reducing noise in the 
results. Since we do not know which background model to 
subtract, we may take as a reasonable guess that the data 
should be close to ACDM and use d L (z) 9 = d L (z) AC ° M as a 
first approximation and then use a boot-strapping method 
to find successively better guess models. Having obtained 
the smoothed luminosity distance, we differentiate it once 
to obtain the Hubble parameter, 
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dz 
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and once again to obtain the equation of state of dark energy 
w( z). 



o(z) 
w(z) = 



[2(1 + 4/3] H'/H - 1 
- {H /H)*n 0m (1 + z 



(3) 



In any kind of smoothing scheme for the luminosity 
distance, some bias is introduced both in d^ and in de- 
rived quantitie s like H(z) and w(z) (see appendix Al in 
IShafieloo etahl l|2006h to find detailed calculations of the 
bias). It is important to choose a value of A which gives 
a small value of the bias and also reasonably small errors 
on derived cosmological parameters. To estimate the value 
of A in {TJ, we consider the following relation between the 
reconstructed results, quality and quantity of the data and 
the smoothing parameters. One can show t hat the relative 
error bars on H(z) scale as (|Tegmarkl l2002) 



5H_ 
H 



jVi/2A 3 / 2 
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where N is the total number of supernovae (for approxi- 
mately uniform distribution of supernovae over the redshift 
range) and a is the noise of the data. From the above equa- 
tion we see that a larger number of supernovae or larger 
width of smoothing, A, will decrease the error bars on 
the reconstructed H, but as it has been reported earlier 
ijShafieloo et al.l l2006h . the bias of the method is approxi- 
mately related to A 2 . This implies that by increasing A we 
will also increase the bias of the results. If we attempt to es- 
timate A such that ^S- oc 3<j, then for N = 182 data points 
(which is the number of data points in the Gold sample), 
we get A = 0.084 for a single iteration of our method. How- 
ever, with each iteration, the errors on the parameters will 
increase. Therefore, using this value of A when we use an 
iterative process to find the guess model will result in such 
large errors on the cosmological parameters as to render the 
reconstruction exercis e meaningless. It has been shown in 
IShafieloo et alj l|2006l ) that at the M-th iteration, the error 
on lnclL will be approximately ^(lndt) ~ \i r M5o(\ndL), 
n\ and the error on lndi, scales as 1/A. We would like the 
errors after M iterations to be commensurate with the opti- 
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mum errors obtained for a single iteration, So, so we require 
A pti m ai — VMAo. Therefore, if we wish to stop the boot- 
strapping after 50 iterations, then A opt i ma i ~ 0.6. However, 
after this rough estimation of the values of A and M, we can 
still play around these values to find the best combination, 
by minimising the likelihood of the reconstructed results to 
the data. In the following, we use A = 0.6 and we calculate 
the x 2 of the reconstructed distance moduli to the data af- 
ter each iteration, and we stop the boot-strapping process 
after reaching the minimum value of x 2 - This effect, that \ 2 
of the reconstructed results goes to a minimum value and 
increases again with iteration is a reflection of the problem 
of some iterative reconstruction algorithms which are not 
error-sensitive. In these cases the noise will be added to the 
reconstructed results after certain number of iterations and 
the iterative process should be stopped after reaching the 
minimum value of \ 2 to get the best result. Similar effect has 
been reported and studied in the Richardson-Lucy deconvo- 
lution algorithm to reconstr uct the form of primordial po wer 
spectrum from CMB data (|Shafieloo fc Souradeerj|l2004f ). 

In appendix A we show that the results are not sensitive 
to the chosen value of A and also to the assumed initial guess 
model. 



3 RESULTS FROM THE GOLD DATASET 

The recently released Gold sample (|Riess et al.l l2006). con- 
sist of 182 supernovae type la which have been gathered 
from five different subsets of data, observed during the last 
16 years. The range of redshift for these supernovae are be- 
tween 0.024 and 1.75. In this section we use this dataset to 
reconstruct h(z), estimate the value of Ho m , and then recon- 
struct w(z). We choose a flat ACDM model with Qom = 0.30 
as the initial guess model in our calculation and we fix the 
value of A (width of smoothing) to be 0.6. After each iter- 
ation, we compute the \ 2 an d we stop the boot-strapping 
process once % 2 reaches its minimum value. 

The x 2 a t an y iteration is calculated from the formula, 

Xrec,j\n-0) — 2jj -g (t>) 

and is followed by marginalising over Ho- We have 
marginalised over Ho by integrating over the probability 
density p oc exp(—x 2 /2) for all values of Ho- 

In Eq[5] /Vec,j(-Ho, Zj) is the reconstructed result at the 
jth iteration for the distance moduli at redshift Zi, assuming 
th e value of Hp , and Mobsfe) is the Gold sample data given 
bv iRiess et al] (2006). In Figure [1] we show the x 2 °f the 
reconstructed results at different iterations, after marginal- 
ising over Ho- As we see, the x 2 has a minimum around 
j — 89 and after this, x 2 ls slowly increasing. So we stop the 
boot-strapping process at this iteration and determine h(z). 
We can also see that for the initial guess ACDM model, the 
Ax 2 of the best recovered result is less than 4 which means 
that the flat ACDM model is in agreement with the Gold 
sample within 2a. 

By marginalizing over the Hubble parameter, we carry 
out a similar trea tment on the data as it has been done by 
IRiess et "all (|2006l ) to calculate the x 2 f° r different cosmolog- 
ical models. As the Gold data are based on a Hubble param- 
eter of 65 km/sec/Mpc, the reconstruction method should 
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Figure 1. Computed \ 2 f° r the reconstructed results at each 
iteration, using Gold sample. 
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Figure 2. Probability density of the best reconstructed result 
from Gold data for different values of Hubble parameter. 

be able to recover this value for the Hubble constant. In 
fact the peak of the probability density of the reconstructed 
result for different values of the Hubble parameter should 
be close to Ho = 65 km/sec/Mpc. In Figure [2] we show 
the probability density of the best reconstructed result from 
Gold data for different values of the Hubble parameter. We 
see that the probability density has a sharp peak around 
Ho =65 km/sec/Mpc. 

We should also note here, that the reduced x 2 of the 
reconstructed results seems to be consistently below 1 (how- 
ever it is not trivial to define the degree of freedom in our 
smoothing method and hence the reduced x 2 , but we can see 
that the resultant x 2 of the reconstructed results is around 
25 less than the number of data points). We can also see in 
Figure [TJthat the reduced x 2 of the first initial guess model, 
which is a ACDM model, is also below 1. It shows that the 
error-bars of the supernovae data points are quite large and 
many different reconstructed results may have a reduced x 2 
of less than 1. In this paper we only calculate the x 2 of the 
reconstructed results and we compare different results by 
calculating the Ax 2 to the best result with a minimum x 2 ■ 

In Figure [3] left panel, we show the reconstructed h(z) 
for the Gold data set. The red solid line has the highest 
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Figure 3. Reconstructed h(z) (left) and q(z) (right) by using Gold dataset. Red solid line is the best recovered result and the green 
dashed lines are within 1<t away from the best result. Based on our results, the transition between deceleration and acceleration phases of 
the universe occurs at 0.38 < z < 0.48 within lcr error-bar from the best recovered result. I n the left panel we can also see 3 uncorrelated 
and independent measurements of h(z) from the Gold sample (blue dotted crosses from iRiess et all (|2006h l for comparision with our 
reconstructed results. 
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Figure 4. Left panel: The derived value of A/\/£lo m from supernovae Gold data within its lcr error-bars (red solid line and green 
dashed lines) in comparison with its measured value from observation of LRGs within its lcr error-bars (blue dotted lines) for different 
values of Qo m . Right panel: Reconstructed w(z) for the Gold dataset. Red solid line is the best recovered result and the green dashed 
lines are within lcr away from the best result. To get these results, we have marginalised over Qom = 0.277 ± 0.022. 



likelihood and is our best reconstruction. All the other lines 
are within la away from the best recovered result. These 
lines are recovered results from our smoothing method by 
using different numbers of iterations in the boot-strapping 
process. The Ax 2 for all of these lines is less than 1, and so 
we can consider them to lie within lcr of the best result. We 
should note that these green dashed lines in Figure [3] are 
in fact a non-exhaustive sample of results which are within 
lcr away from the best recovered result. As we see in the 
Figure [3] the reconstructed h(z) at high redshift has a very 



big degeneracy. This is expected since there is only a single 
supernova beyond redshift 1.4! 

In this figure we can also see 3 uncorrelated and inde- 
pendent measurem ents of h(z) from t he Gold sample (blue 
dotted crosses from lRiess et all l|2006h ) for comparison with 
our results. We can see that these two results are consistent 
with each other within t heir la limits. However w e should 
mention here that in the jWang fc Tegmarkl (2005) method 
used bv lRiess et al.l (|2006h for uncorrelated estimates of the 
expansion history, there is a slight bias in the reconstruc- 
tion of h(z) and the higher derivatives of the data. It is 
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mainly because of using the average of the measured quan- 
tity ^(z)" 1 , which typi cally is not a straight lin e. So as it 
has been mentioned in IWang fc Teemarkl 1 2005 ) . the mea- 
sured average of h(z)~ l (and hence h(z)) over a redshift bin 
will generally lie either slightly above or below the actual 
curve at the bin center. This can be the reason that why 
the centres of the crosses in Figure [3] left panel, for uncorre- 
cted estimates of the expansion history, are slightly above 
or below our reconstructed curve for the h(z). 

To reconstruct the Hubble parameter, h(z), we do not 
need to know the value of Qom- Another important cosmo- 
logical quantity which we can derive from the reconstructed 
h(z) (independent of the value of Qom), is the deceleration 
parameter, q(z), 



q(z) = (l + z) 



H'( 



H(z) 



(6) 



In Figure [3] right panel, we show reconstructed q(z). 
For the Gold data our method shows that the transi- 
tion between deceleration and acceleration occurred at 
0.38 < z < 0.48 (at la). The best reconstruction shows 
the redshift of transition to be z a — 0.42. This is in 
agreement with results obtained using parametric meth ods 
jAlam. Sahni fc Starobinskvll2007l ; iGong fc Wan3l200rj ). 

To derive the equation of state of dark energy w(z), 
one needs to know the value of fiom, as we see in EqO To 
estimate the value of £7o m , without using any parameteri- 
sation and in a model-independent way, we can use the re- 
sults of the detection of the baryon a coustic oscillation peak 
(|Linde r 2003; Eiscns tein et al.ll2005l ). The distance factor A 
up to redshift 0.35, measured by observation of luminous 
red galaxies in detection of baryon acoustic oscillation peak 
(which have been claimed to be relatively independent of the 
model of dark energy), can be derived directly for different 
values of Qom by using the reconstructed h(z), 



A 



vHom 



i 



Z1 dz_ 
h{z) 



2/3 
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where the measured value of A is A — 0.469( ^ )-°- 35 ± 
0.017 at 2i = 0.35. The 3-year WMAP results, when com- 
bined with the results of baryon acoustic oscillations, yield 
n = 0.951 for the spectral index of the primordi al power 
spectrum l|Spergel et al.ll2006l ; lLAMBDA web site! '). By us- 
ing the best reconstructed results for h(z), we get A/y/fl m = 
0.901. In Figure [4] left panel, we see the derived value of 
A/\/Q. m from supernovae data in comparison with its mea- 
sured value from observation of LRGs for different values 
of fiom- It is clear that these two independent observations 
which are completely different by nature, are very much in 
agreement if 0.255 < fiom < 0.299. This derived value of 
flom is completely independent of any dark energy model as- 
sumption (within the framework of standard general relativ- 
ity) and is in very close agreement with t he results from large 
scale structure measurements from 2dF (|Tegmarkll2o"ol ) and 
SDSS (jCole et al.ll2005h . This derive d value of Qp m is also i n 
good agreement with the results from lFav fc Tavakoll (2006), 
where a different model-independent method of reconstruc- 
tion has been used. 

Now by marginalising over Qo m = 0.277 ± 0.022, which 
is the range of agreement between the two observations, 
we can reconstruct w(z) from our previously reconstructed 
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Figure 5. Computed x 2 f° r the reconstructed results at each 
iteration for the SNLS datasct. 



h(z). In Figure [4] right panel, we show the reconstructed 
w(z), marginalised over Oo m for the Gold dataset. We see 
that the data prefer evolving dark energy to the cosmolog- 
ical constant. The degeneracy for the equation of state of 
dark energy at high redshifts is very large and it is almost 
impossible to say much about w(z) at high redshifts. 



4 RESULTS FROM THE SNLS DATASET 

In this section we use the same procedure as we used in 
the previous section to deal with SNLS supernovae data. 
The SNLS dataset contains 115 data points in the range of 
0.1 < z < 1.0. We use this dataset, first to reconstruct the 
Hubble parameter, h(z), and the deceleration factor, q(z), 
up to redshift 1. Then by using the results of detection of 
baryon acoustic oscillation peak we derive the value of fiom , 
following which we recover the form of w(z). We use the 
distance modules of the su pernovae available in Tables 8 
and 9 in lAstier et all (2005) as our dataset in this section. 

In Figure \5\ we see the computed \ 2 f°r t ne recon- 
structed results using smoothing method at each iteration. 
As we see, the \ 2 diverges to its minimum value very fast at 
just the 5th iteration. In Figure[6]we show the reconstructed 
h(z) (left panel) and q(z) (right panel) for the SNLS dataset. 
The red solid line has the best likelihood, which is our best 
reconstructed result. All the other lines are within la away 
from the best recovered result. We should like to emphasise 
here that these results (green dashed lines) are not represen- 
tative of all the possibilities which give the likelihood within 
la of the best recovered result. However they can show the 
overall behaviour of the quantities which we have studied. 
Our results for SNLS data show that the transition from 
deceleration to acceleration phase of the universe occurs at 
redshifts higher than 0.7. The fact that we cannot put an 
upper limit to the redshift of the commence of acceleration 
is due to the absence of supernovae data at z > 1 in SNLS 
dataset. 

As we have discussed earlier in the previous section, we 
use the results of detection of baryon acoustic oscillation 
peak to determine the value of Qom ■ Then by marginalising 
over the recovered value of flo-m, we derive the dynamics of 
w(z). In Figure [7] we see the derived value of Qom and the 
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reconstructed form of w(z). We see that the ACDM model is 
in much better agreement with SNLS data than with Gold 
data. 

By comparing the recovered results from the SNLS and 
Gold datasets, we can clearly see an inconsistency between 
these two supernovae datasets. This inconsistency is obvi- 
ous by looking at the reconstructed q(z) and w(z) in the 
middle and high redshift ranges. Gold data suggest the red- 
shift of the commencement of acceleration at z a ~ 0.42 
while SNLS data suggest z a ~ 0.80. The reconstructed w(z) 
from these two datasets also shows a very different behav- 
ior in the middle and high redshift ranges. The discord be- 
tween Gold and SNLS super novae datasets has been re- 
ported and studied earlier by iNesseris fc Perivolaropoulosl 
(2006), and a similar discord between Gold supernovae 
data and other cosmological observations like Cosmic Mi- 
crowave Background observations from WMAP and observa- 
tions of cluster abundance, has a lso be en reported earlier by 
IJassal. Bagla fe Padmanabhanl (|2006T ). However SNLS su- 
pernovae data seem to be in good agreement with the other 
cosmological observations. Based on all these results and 
analyses, we may conclude that some significant systemat- 
ics in the Gold data (or in a part of the data) might be the 
reasos for these inconsistencies. 

Interestingly, the recovered values of Qom from Gold 
and SNLS data (by using the results of detection of baryon 
acoustic oscillation peak), are in very close agreement. In 
both cases the derived value of Qom is around 0.276 ± 0.022. 
We should note here that the two data sets rely on pretty 
much the same nearby supernovae samples and that is why 
the results are similar in this range. It is something which 
we logically expect to get. But in fact it shows one of the ad- 
vantages of this method over the functional fitting methods. 
By using a functional fitting method, the recovered results 
at any redshift would be equally dependent on the data set 
in the whole redshift range. But here, by using our smooth- 
ing method, we can clearly see that despite the significant 
differences between the reconstructed results from Gold and 
SNLS datasets in the middle and high redshift regions, the 
reconstructed results for the expansion history at low red- 
shifts (which we use to estimate the value of matter den- 
sity), are not affected by the big differences between the two 
datasets at the higher redshifts. 



and SNLS datasets. This determination is directly related to 
the supernovae data points at redshifts lower than z = 0.35. 
We have got Q 0m ~ 0.276 ± 0.023 for both Gold and SNLS 
datasets, which is in good agreement with results of SDSS 
and 2dF large scale structure observations, and also with 
results of recent Chandra X-ray o bservations of the relaxed 
galaxy clusters (|Allen et al.ll2007l ). 

This derived value of Qo m also agree with the recent 
WMAP 3 years CMB data, if we assume the broken scale 
invariant spectrum for the fo rm of the primordial spectrum 
l|Shafieloo fc Sourade cp 2007). In the derivation of q(z) and 
the stage of transition from deceleration to acceleration in 
the dynamics of the universe, we found disagreement be- 
tween Gold and SNLS datasets. Gold data suggest the red- 
shift of the commence of acceleration at z a — 0.42 while 
SNLS data suggest z a ~ 0.80. 

After marginalizing over the derived value of Qom, we 
have reconstructed w(z). The inconsistency between Gold 
and SNLS supernovae datasets is also obvious by looking at 
the reconstructed w(z) from these two datasets. The derived 
form of w(z) from SNLS dataset, is in good concordance 
with ACDM model, while Gold dataset prefers an evolv- 
ing form of dark energy (however ACDM is still in agree- 
ment with the Gold dataset to within 2a). This discrep- 
ancy between Gold an d SNLS datasets has been r e porte d 
earlier by other groups ((Nesseris fc Perivolaropoulosl {2006); 
lAlam. Sahni fc Starobinskvl ((2007)). As the Gold sample is 
also relatively in disagreement with the other cosmological 
observ ations like CMB and observations of cluster abun- 
dance l|Jassal. Bagla fc Padmanab han (2006)), we may con- 
clude that the effect of systematics in the Gold dataset (or 
at least in a part of the data) is significant. 

The large error-bars at the high redshifts for the recon- 
structed results, reflect the significant lack of data points. 
This effect may not be seen if we use some of the parametric 
methods of analysis, but as we deal with the data directly 
here, we notice that the lack of data points at high redshifts 
limits our ability to say much about the behavior of the 
Universe at the early stages at high redshifts. This is an- 
other important feature of our smoothing method in which 
the reconstructed results at any redshift rely mostly on the 
supernovae data points at the same redshift range. 



5 DISCUSSION AND CONCLUSION 

In this paper we have shown t hat by improving the effi- 
cacy of the smoothing method l)Shafieloo et al.l 20061 ) . we 
can reconstruct the expansion history of the universe in 
a model-independent way, using current supernovae data. 
We have used the smoothing method to reconstruct the ex- 
pansion history of the universe, h(z), the deceleration pa- 
rameter, q(z), the value of flo m and the equation of state 
of dark energy, w(z), independently of any assumption of 
the theoretical model of the universe, within the framework 
of standard general relativity. This is an advantage of this 
method over the functional fitting methods where the re- 
sults are usually biased by the form of the functional fitting 
or the assumed theoretical model. We dealt with two recent 
datasets, Gold and SNLS in our analysis. In determining the 
value of Oo m , we found excellent agreement between Gold 
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APPENDIX A: EXAMINING THE 
ROBUSTNESS OF THE METHOD 

In this section we show that the results of the smoothing 
method are robust against the choice of the initial guess 
model and also to the chosen value of A. 

We assumed three different cosmological models as our 
initial guess model and we applied our smoothing method on 
the Gold dataset. The final results by using these three dif- 
ferent initial guess models are almost identical with A^ 2 < 
0.01. We have got \ 2 = 157.40 by using a flat ACDM model 
with Qom ~ 0.30 as the initial guess model after 89 iteration, 
while we have got x 2 = 157.40 for a flat ACDM model with 
Ho™ = 0.25 after 91 iteration, and \ 2 ~ 157.39 for a flat 
quiessence model with Qom = 0.30 and w(z) = —0.8 after 
104 iteration. In figure [Ail we can see the reconstructed h(z) 
and w(z) for the Gold data set by assuming these three dif- 
ferent initial guess models. As we see, the robustness of the 
method for the choice of the initial guess model is obvious. 

We have also used different values of A (width of 
smoothing in Eq[l}, in our reconstruction process to check 
the reliability and stability of our results against the changes 
in the value of A. We have used three values of A equal to 
0.30, 0.60 and 0.90 in our smoothing method and we have ap- 
plied it on the Gold dataset. By using A = 0.30 we have got 
X 2 = 157.38 after 9 iteration, while we have got x 2 = 157.40 
by using A = 0.60 after 89 iteration, and x 2 ~ 157.41 by 
using A = 0.90 after 407 iteration. In figure lA2l we can see 
the reconstructed h(z) and w(z) for the Gold dataset by 
using these three values of A. We can clearly see that the 
results are not sensitive to the given value of A. These two 
examinations confirm the overall robustness of the method 
for different initial assumptions. 
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Figure Al. Reconstructed h(z) (left panel) and w(z) (right panel) for the Gold dataset by assuming three different initial guess models. 
The red solid line is the reconstructed result by using a flat ACDM model with Qom = 0.30 as the initial guess model. The green dashed 
line is the reconstructed results by using a flat ACDM model with Qom = 0.25, and the blue dotted line is the reconstructed result by 
using a fiat quiessence model with w(z) = —0.8 and Qo m = 0.30 as the initial guess models. We can clearly see that the results are 
almost identical which shows the robustness of the method for the different choices of the initial guess model. 
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Figure A2. Reconstructed h(z) (left panel) and w(z) (right panel) for the Gold dataset by using three different values of A (width of 
smoothing). The red solid line is the reconstructed result by using A = 0.60. The green dashed line is the reconstructed results by using 
A = 0.90, and the blue dotted line is the reconstructed result by using A = 0.30. In all these cases we have stopped the boot-strapping 
process after reaching to the minimum \ 2 ■ We can see that the method is robust against the variation of A in a wide range. 



